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O 

The quark behaviour in the background of intensive stochastic gluon field is studied. An approximate procedure 
for calculating the effective Hamiltonian is developed and the corresponding ground state within the Hartree- 
Fock-Bogolyubov approach is found. The comparative analysis of various model Hamiltonian is given and 
transition to the chiral limit in the Keldysh model is discused in detail. We study response to the process of 
filling up the Fermi sphere with quarks, calculate the vacuum pressure and demonstrate the existence of filled- in 
state degenerate with the vacuum one. 

' We study the quark (anti-quark) behaviour while being influenced by intensive stochastic gluon field and 

work in the context of the Euclidean field theory. The corresponding Lagrangian density is the following 

Oh" 

CD ' Ce — q {iluDu, + im) q , (1) 

here q (q) — are the quark (anti-quarks) fields with covariant derivative = — igA^t a where A°\ is the 
gluon field, t a = A a /2 are the generators of colour gauge group SU(N C ) and m is the current quark mass. 
■ As the model of stochastic gluon field we refer to the example of (anti-)instantons considering an ensemble 

of these quasi-classical configurations. On the way to construct an effective theory we consider the quenched 
approximation and neglect all the contributions coming from gluon fields A ex generated by the (anti-)quarks 
A ex -C A. Then the corresponding Hamiltonian description results from 

n = nq-C E , n = ^=zq+, (2) 

o\ : _ q 

and H,q = —q (17V + im)q, for noninteracting quarks. In Schrodinger representation the quark field evolution 
is determined by the equation for the quark probability amplitude ^ as 



1 — 1 

X 



* = -H^ , (3) 

with the density of interaction Hamiltonian Vs = <z(x) i°7 M A^(i,x) g(x). The explicit dependence on "time" is 
present at the gluon field only. The creation and annihilation operators of quarks and anti-quarks a + , a, b + , b 
have no "time" dependence and consequently 

<?m(x) = J (2\ P l\) 172 [ a (P' s ' c ) u ai (p,s,c) e !px + 6+(p,s,c) v ai (p, s,c) e _4px ] . (4) 

The stochastic character of gluon field (which we supposed) allows us to develop the approximate description 
of the state "J if the following procedure of averaging ^ — > (\f) = J * dr ^(r)/ t is intoduced. With this 

procedure taken the futher step is to turn to the approach of constructing a density matrix However, 
here we believe that at calculating the ground state (or more generally with quasi-stationary state) it might be 
sufficiently informative to operate with the averaged amplitude directly. Then in the interaction representation 
* = e ff °*$ we have the equation for state $ as $ = — V<&, V = e Hot Vse~ Hat . Now the "time" dependence 
appears in quark operators as well and after averaging over the short-wavelength component one may obtain 
the following equation 

/•OO 

($(«)) = +/ dr (V(t)V(t-r)) (<Z>(t)) . (5) 







The limitations to have such a factorization validated are well known in the theory of stochastic differential 
equations (see, for example, PQ). The integration interval in Eq.Q may be extended to the infinite "time" 
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because of the (supposed) rapid decrease of the corresponding correlation function. Now we are allowed to 
deal with amplitude (3>(f)) in the right hand side of Eq.f5|) instead the amplitude with the shifted arguments 
in order to get an ordinary integro-differential equation. In the quantum field theory applications it is usually 
difficult to construct the correlation function in the most general form. However, if we are going to limit our 
interest by describing the long-wavelength quark component only then gluon field correlator (A a ^{x)A b v {y)) may 
be factorized and as a result we have 



($(*)) = jdxg(x,t) f 7fl g(x,t) J dr J dy q(y,t-r) t b lv q(y,t-r) ^(i;(i,x)^(i-r,y)) <$(*)) . 

Having assumed the correlation function rapidly decreasing in "time" we could ignore all the retarding effects 
in the quark operators. Turning back to the Schrodinger representation we have for the state amplitude x — 
e~ Hot (Q) the following equation 

X = -H lndX -, H ind = -q (*7V + im) q - q t a 1{1 q J dy q' t b lv q' dr g 2 {A^A b ) , (6) 

with q = q(x), q = g(x), q' = q(y), q' = q{y) and At = A°(t,x), A b — A h v {t — r, y). Now the correlation 
function might be presented as J °° dr g 2 {A®A b ) = S ab 8^ J(x — y) + J M „(x — y). In our consideration we 
ignore the contribution of the second formfactor spanning on the components of the vector x — y. Thus, on 
output we receive the Hamiltonian of four-fermion interaction with the formfactor rooted in the presence of 
two quark currents in the points x and y. With this form of the effective Hamiltonian we could apply the 
Hartree-Fock-Bogolyubov method to find its ground state as one constructed by the quark-anti-quark pairs 
with the oppositely directed momenta 

\a)=T\0), T = n p , SiC exp j °- [ o+(p, s, c) 6 + (-p, s, c) + o(p, s, c) 6(-p, s, c) ] J , (7) 

where the parameter 9(jp) characterizes the pairing strength. Minimizing the mean energy functional one is able 
to determine the angle 9 mag nitude i<£i^iM = o. By the help of dr essing tranformation T we introduce the 
creation and annihilation operators of quasi-particles A ~ T a T -1 , B = T b T _1 . Dropping the calculation 
details out (see Ref. [2}) we present here the following result for the mean energy as a function of the 9 angle 

/ i it i \ [ dp 2N c pi 

_G (-(37- J)J^L + (47- J)-^f - ( S m6-™ cose) (sin9' - ^cos 

J (2tt)6 \ > 4 ||<? 4 | 1^4 1 1^ 4 | V P J V 9 

ui \ I m 



Pi gj 

\pi\\n\ V " ' p "J V" ' i 



-(-2I5 i:j - 2Jij + JSij) f- J j ( cos6» + — sin6» ) ( cos^' + — sin( 



here the following designations are used p = |p|, q = |q|, / = I(p + q), Jij = Jij(p + q), J = x Ju, 
p 2 = q 2 = —rn 2 , 9' = 9{q) where G is the constant of corresponding four-fermion interaction (the relevant 
details can be found in [2]). The first integral in Eq. ([5]) comes from free Hamiltonian, and we make a natural 
subtraction (adding the unit) in order to have zero mean free energy when the angle of pairing is trivial. 
Nambu— Jona-Lasinio model. In order to get an idea of the parameter scales we continue with handling 
the model in which the formfactor behaves in the coordinate space as /(x — y) = 5(x — y), J^ v — 0, dropping 
contribution spanned on the piqj tensor also. Actually, it corresponds to the Nambu-Jona-Lasinio model [3]. 
As well known the model with such a formfactor requires the regularization and, hence, the cutoff parameter A 
comes to the play 



W 



A dp 



(2tt) 3 



Ip^l (1 — cos 9) — G-r^— ( sin 6* cos 0] f ^ ( sin 9' cos( 

m\ \ P J J (2tt) j \q 4 \ V q 



0) 



We adjust the NJL model with the parameter set given by Hatsuda and Kunihiro [3] in which A = 631MeV, 
to = 5.5MeV. One curious point of this model is that the solution for optimal angle 9 in the whole interval 
p G [0,A] can be found by solving the simple trigonometrical equation (p 2 +m 2 ) sin (9 — M q (pcos9 + msm9) = 0, 

with the dynamical quark mass M q = 2G J A ( 2 tT) 3 JpT\ ( sm ^ — ^ cos 6*^ . Eventually the results obtained look 
like M q = —335 MeV for dynamical quark mass and (a\qq\a) = — i (245 MeV) 3 for the quark condensate with 
the following definition of the quark condensate {a\qq\a) = J °° dp {psva9 — to cos 9). The Keldysh 
model. Now we are going to analyse the limit, in some extent, opposite to the NJL model, i.e. we are dealing 
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Figure 1. Phase portrait of the Keldysh model, 
sin 9 as a function of momentum p(MeV) (dashed 
line for imaginary values). The dotted curves cor- 
responds to the solution in the chiral limit rn = 0. 



Figure 2. The optimal angle 9 as a function of 
momentum p(MeV) . The solid line corresponds to 
the NJL model and the dashed one to the Keldysh 
model. The current quark mass is m = 5.5 MeV 
and p e ~ 40 MeV. 



with the formfactor behaving as a delta function but in the momentum space (analogously the Keldysh model, 
well known in the physics of condensed matter [4]), I(p) = (2tt) 3 <5(p). Here the mean energy functional has 
the following form 



W(m) = 



dp 

(2VF 



fP" ( Tfl 

\pi\ (l — cos8) — G- — — \sm9 cos( 

m\ V P 



(10) 



contrary to the NJL model there is no need to introduce any cut off. The equation for calculating the optimal 
angle 9 becomes the transcendental one \p4\ 3 sin 9 — 2G (pcos8 + m sin 9) (psin^ — m cos 9) = 0, and, clearly, 
it is rather difficult to get its solution in a general form. Fortunately, it is much easier and quite informative to 
analyse the model in the chiral limit m = 0. There exist one trivial solution 9 = and two nontrivial ones (for 
the positive and negative angles) which obey the equation cos 9 = ^ . Obviously, these solutions are reasonable 

if the momentum is limited by p < 2G. Then for the mean energy (for real solution) we have: W±(0) = —j§^2, 

and for the quark condensate: {a\qq\cr)(0) = 1 ■ For the trivial solution the mean energy equals to 



zero together with the quark condensate Wo(0) = 0, (o-\qq\a) (0) = 0. Introducing the practical designation 
sin# = ^ p2+ j )^2)i/2 which characterizes the pairing strength by the parameter Mg we have, for example, for the 

1/2 

nontrivial solution Mg = (AG 2 — p 2 ) . In order to compare the results with the NJL model we fixed the value 
of four-fermion interaction constant as Mg (0) = 2G = 335 MeV. It is interesting to notice that the respective 
energy becomes constant E(p) = ^Jp 2 + M| , E(p) = 2G. 

After having done the analysis in the chiral limit which is shown by the dotted line in Fig. 1 we would like to 
comment the situation beyond this limit. The evolution of corresponding branches is available on the same plot 
[TJ The minimum of mean energy functional can be realized with the piecewise continuous functions. At the 
local vicinity of coordinate origin we start with some branch of the solution, then relevant solution passes from 
one possible branch to another one at any subinterval. But in any case there is only one way to continue the real 
solution at streaming to the infinite limit. As to the functional (TIT))) the contribution of the term proportional 
to the cosine in the second parenthesis is divergent even if the angle 9 is zero. It means the mean energy 
out of chiral limit goes to an infinity at any nonzero value of quark mass. The same conclusion is valid for 
the chiral condensate. In principle this functional could be regularized and corresponding continuation might 
be done but it is out of this presentation scope (see Ref. [5]). It is not difficult to demonstrate the similar 
discontinuities of functional are present, for example, for Gaussian 7(x) = G exp (— a 2 x 2 ), and exponential 
I(x) = G exp (—a |x|), formfactors and they are present even in the NJL model but this fact is masked by the 
cut off parameter. Comparing the optimal angles in the NJL and Keldysh models (see Fig. 2) it is interesting 
to notice that the formation of quasiparticles becomes significant at some momentum value close to the origin 
pg ~ 40 MeV (for the Gaussian and exponentional formfactors it is around pg ~ 150 MeV) but not directly at 
the zero value. It is clear the inverse value of this parameter determines the characteristic size of quasiparticle. 
Analysing the discontinuity of mean energy functional and quark condensate we face some troubles at fitting the 
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Figure 3. The quark chemical potential as 
a function of the Fermi momentum for the NJL 
model. The solid line corresponds to the current 
quark mass m = 5.5 MeV and the dashed one 
shows the behaviour in the chiral limit. 



Figure 4. The pressure of the quark ensemble 
as function of the Fermi momentum. The solid 
line for the current quark mass m = 5.5 MeV. The 
dashed one in chiral limit. 



quark condensate, for example. However, the dynamical quark mass and quark condensate are nonobservable 
quantities and it is curious to remark here that although the mean energy of the quark system is minus infinity 
the meson observables are finite and even in Keldysh model the mesons are recognizable with reasonable scale 
and we can in principle make a fit for this observables [5] . 

Now our central issue could be formulated in the following way — to construct the state filled in by quasi- 
particles (the Sletter determinant) \N) — nipi<Pp-g A + (V\s) \<j), which possesses the minimal mean energy 
(N\H\N) (surely, we assume the quasi-particles are stable). Here Pp stands for the Fermi momentum and 
the polarization runs over all possible values. It allows us to optimize the dressing transformation and, as the 
consequence, to follow up the modifications of quasiparticles being influenced by the process of filling in the 
Fermi sphere. Eventualy it fixes the form of charge operator (particle number operator) |(7V|gi7 4 (7|iV}|. Let us 
define the partial energy density per one quark degree of freedom, as w — £ = E/V where E is the total 
energy of ensemble. For the ensemble of quasi-particles we obtain the following expression for the partial energy 

- G / 7^3 i — i sm 9 cos9 l / 7o~TT , — r sm cos J . 

Jp f (2it) 3 \p 4 \ V P /7p f (2tt) 3 \q 4 \ V Q J 

It could have ruther interesting interpretation if compared to the vacuum mean energy Eq.®. It is easy to 
see that for the state with the filled-in Fermi sphere the angles of pairing could be defined by the condition of 
functional minimum only for the momenta larger than Fermi momentum Pp . Then the quarks composing 
the Fermi sphere look like the free (non-interacting) ones, as seen from the first term of Eq. (fTTj) . Now let 
us calculate the quark chemical potential which, by definition, is an energy necessary for adding (removing) 
one quasi-particle to (from) a system \i = where N — 2N C V J F j 2 R$ — V Pp is the total number 
of particles in the volume V. Redefining the chemical potential as ji = j?r -gp^ we consider the model with 
correlation function behaving as the (5-function in the coordinate space, which corresponds to NJL model. The 
following relation could be obtained in this case (see Ref. [6j) 

M =[P| + M 9 2 ] 1 / 2 . 

Let us remind that for the free fermion gas the chemical potential increases monotonically with the Fermi 
momentum growing. The curious feature of the NJL model is the appearance of state almost degenerate with 
the vacuum state while the process of filling up the Fermi sphere reaches to the momenta close to the dynamical 
quark mass value (the similar value is peculiar to the momentum of quark inside a baryon), see Fig. [3l This 
state density with the factor 3 (which expresses the relation between baryonic and quark degrees of freedom) 
absorbed corresponds to a normal nuclear density (n ~ 0.12/fm 3 ), and chiral condensate could be estimated 
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as | < qq > 1 / 3 | ~ 100 MeV. In the chiral limit the chemical potential is close to the discussed point and 
is even smaller than the vacuum one. The full coincidence of the chemical potentials occurs at the values of 
current quark mass around 2 MeV. In fact, Fig. [3] shows that the u quark bond looks stronger than one of the 
d quark. The pressure of the quark ensemble P — — ||p = — ^ + ^ = — £ + fj, n, is depicted in Fig. 2] 
as a function of the Fermi momentum where n = N/V is the quark density. The quark pressure at the values 
of the Fermi momentum close to the quantity of dynamical quark mass is approximately degenerate with the 
vacuum pressure (slightly lower than the vacuum one). The vacuum density is of order 40 — 50 MeV/fm 3 and 
corresponds well to the value extracted from the bag models. Apparently our estimate of the effects responding 
to the process of filling up the Fermi sphere entails a hope to understand a routine feature of hadron world, 
namely, the fact of quark equilibrium in the vacuum and inside the proton. The chemical potential degeneracy 
and specific behaviour of the quark pressure (with one new essential element which is just the presence of 
instability region dP/dPp < 0) justify, in principle, the conventional bag model. It urges to consider the filled 
states \N) as natural 'building' material for baryon octet (on the strong interaction scale only). 
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